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ABSTRACT 



Aims. We present synthetic spectra around maximum for the type II supernova SN 2003Z, which was first detected on January 29.7 
2003. Comparison with observed spectra aim at the determination of physical parameters for SN 2003Z. 

Methods. Synthetic spectra are calculated with our stellar atmosphere code PHOENIX. It solves the special relativistic equation of ra- 
diative transfer, including large NLTE-calculations and line blanketing by design, in 1 -dimensional spherical symmetry. The observed 
spectra were obtained at the 3.5 meter telescope at Calar Alto. The TWIN instrument was used so that a spectral range from about 
3600 to 7500 A was covered. The spectra were taken on Feb. 4, 5, 9, and II, 2003. 

Results. The physical parameters of the models give the luminosities, a range of possible velocity profiles for the SN, an estimate of 
the colour excess, and the observed metalicity. 
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1. Introduction 

Type II supernovae (SNe II) are thought to originate from the 
core collapse of massive stars. By definition SNe II show strong 
B aimer lines in their spectra, and thus are thought to be from 
stars with much of their hydrogen envelopes intact at the time 
of the explosion. Largely through the use of HST, progenitors 
have been identified for several SNe II, and they seem to come 
from rather low mass stars 9 > M < 12 Mq ('Maund et al."'2005t 
Maund & Smartt 2005; S martt et all 2003; Smartt et al. 2004; 
[Van Dvk et aI.II2003allblld) . 

Since the spectra of SNe II form in hydrogen dominated at- 
mospheres with relatively simple density structures, they should 
be among the most accurate spectra to model with detailed 
model atmosphere codes. 

Indeed, SNe II can be modeled in detail and reddening, 
primordial metalicity, a r id eve n distances can be determined 
jBaron et al. l l2000ll2003L 120041) . Thus, through detailed spectral 
modeling we can determine important physical parameters that 
can be compared with stellar evolution and explosion modeling. 

Here, we us e the PHOENIX relativistic model atmosphere 
code pac kage ('Hauschildt & Baron' l99^, 'HauschildLetin 
119971; [B^on & Hauschildt 1998; Hauschildt et al. 2001) to sim- 
ulate the SN II atmosphere for SN 2003Z during its optically 
thick (photospheric) phase. 

In the following we first describe the observations of SN 
2003Z, and then describe our modeling process including the as- 



sumptions made to construct a model. After that we present our 
results about luminosities, extinction, velocity profiles. Finally 
the results are summarized in the conclusion. 



2. The supernova SN 2003Z 



* Based on observations collected at the Centro Astronmico Hispano 
Alemn (CAHA) at Calar Alto, operated jointly by the Max-Planck 
Institut fr Astronomic and the Instituto de Astrofsica de Andaluca 
(CSIC) 



SN 2003Z was discovered on Janu ary 29, 20 03 (Bol es et akl 
2003) and is a type II supernova (Mat heson et a l. 2003). 

The SN was observed on the 4th, 5th, 9th, and 11th of 
February 2003 with the 3.5 meter telescope at Calar Alto. The 
TWIN spectrograph was used, covering the spectral range from 
about 3600 to 7400 A. The SN is located near the galaxy NGC 
2742 and is believed to lie in the outer part of a spiral arm of 
the galaxy, typical f or a Type II super nova. The redshift of NGC 
2742 is z = 0.0043 (iFalco et al.lll999h . 

From NEE0 th e foreground colour excess of the galaxy is 
E(B - y) = 0.039 (Schlege l et allll998l) . Information about the 
radial velocity of the SN itself or the colour excess of the host 
galaxy are unavailable. 

The physical structure of the atmosphere changes over time 
due to the expansion and the gradual cooling of the ejected mate- 
rial. Therefore, the spectrum also changes over time. From day 
to day the changes are small; however, the spectra of the first 
and the last day of the observations show large differences as 
displayed in Fig.[T] 

The gradual change in the wavelength shift of the spectral 
lines originates from the changing optical depth of the cooling 

' This research has made use of the NASA/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propulsion Laboratory, 
California Institute of Technology, under contract with the National 
Aeronautics and Space Administration. 
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Fig. 1. Comparison of the spectra of the first and last day of the 
observations. A shift in wavelength towards the red and a de- 
crease of flux in the blue part are visible, indicating the cooling 
and slowing of the material producing the spectra. 



and expanding envelope. As time progresses we see into slower 
and "deeper" layers of the atmosphere so that the lines show 
smaller redshifts in the later spectra. 



3. Modeling 

PHOENIX is a general purpose NLTE stellar atmosphere code 
package. For the calculation of the radiative transfer with 
PHOENIX the structure of the atmosphere is approximated by a 
spherically symmetric shell, thereby reducing the system to one 
spatial dimension. Furthermore, the envelope is assumed to be 
homologously expanding, giving a linearly increasing velocity 
field (Sedov solution) 



v(r) = vo 



^0 



(1) 



with v'o being the velocity and Rq the radius at t = 1 in the radial 
optical dep th grid. The run of th e density is approximated by a 
power law (IZel'dovich & Raizerll 19671) 



(2) 



Due to the fact that the atmosphere of a SN is expanding 
at ^ 10 percent of the speed of lig ht, the special relativistic 
equation (iMihalas & Weibel Mihalaail984 ) of radiative transfer 
in moving media must be solved. In co-moving form written in 
terms of wavelengths we have 



with : a;[ - y 



OS OA 



(3) 



or 



where I,i is the specific intensity and the differentiation of the 
path-length s is the differentiation along the monochromatic 
path of a photon in the co-moving frame. Furthermore, rfx is 
the emissivity and xa the extinction coefficient. The emissiv- 
ity T] depends via scattering on the mean intensity J,\ itself, 
therefore, the problem is not a mere differential equation but 
an integro-differential equation. The radiative transfer problem 
is solved with an operator splitting technique in the co-moving 



frame dHauschildt & BaronI 120041: iHauschildtl |1992|) . For the 
NLTE models, the coupled radiation transport and rate equa- 
tion problem is so lved by a multi-level operator splitting method 
(lHauschildtlll993h . 

Details of the code and the numerical metho ds used in 
PHOENIX can be found in lHauschildt & BaronI (Il999h . 

All models were calculated on a radial grid with 100 layers, 
using a logarithmic optical depth grid (r-grid) in the continuum 
at 5000 A with a range from t = 10^^ - 10""'. 

Other model parameters vary from model to model. For in- 
stance the final model for the first day of the observation had a 
pressure at t = 1 of 4.44 ■ 10"** bar and the corresponding radius 
was 9.62 ■ lO'"* cm. 

Unstable nuclei - most importantly ^^Ni - are created in SN 
outburst and, therefore, the input of radioactive energy plays an 
important role in the formation of the spectrum. The treatment of 
the influence of radioactivity is simplified by the assumption that 
y-ray deposition follows the density profile and the total mass of 
radiative nuclei is a parameter in the modeling. In all models 
we iterated the temperature structure to fulfill the condition of 
energy conservation in the co-moving frame. Furthermore, the 
most important species were treated in complete NLTE. (see 
Table [1]). The decision which ionization stages were included 
in NLTE was made by taking into account the partial pressures 
of those stages. An ion was ignored for the NLTE calculation 
only if the pressure of an ion was for all layers was at least 15 
dex smaller than the pressure of the dominant ion of the species. 

To model SN atmospheres, a few physical parameters must 
be specified. For example, the luminosity, the exponent n of the 
density power law (|2]i, and the velocity field of the atmosphere 
must be specified. These parameters are not known in advance. 
Therefore, several models with different physical parameters had 
to be calculated and the best set of parameters was determined by 
comparison of the observed and synthetic spectra. The first ob- 
servation was modeled to determine the basic parameters of the 
SN and the following observations were reproduced by changing 
the physical parameters appropriately. 



4. Results 

To compare the model spectra with the observation a couple of 
effects have to be taken into account, such as the redshift of the 
SN, and extinction. 

Nothing is known about the intra-galactic extinction in NGC 
2742, hence we only correct for the known foreground extinction 
of NGC 2742 E{B - V) = 0.039. 

We correct the model spectra for the redshift z - 0.0043 
of the host galaxy and for the wavelength shift due to the the 
transition from vacuum to aifl 

First we ran a couple of models varying the luminosity 
(model temperature) and the velocity field to get a first guess for 
those basic parameters. After including the first four ionization 
stages of iron as well as hydrogen and helium in complete NLTE, 
a spectrum was obtained that fits some of the key features of the 
observed spectrum, see Fig. |2] The displayed model spectrum 
has Tmodei - 5600K, a density exponent n = 1 1 and vq = 6300 
km s '(see Eq. [T]i. The model spectrum is blue shifted with re- 
spect to the observation. This is due to difficulties with the con- 
vergence of the models that had lower photospheric velocities. 
This initial result is clearly due to a velocity that was too large 
in the models. This simple preliminary model had few species in 



This correction is very small and we apply it only for completeness. 
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NLTE and it was difficult to converge models with lower veloc- 
ities. As we improved the input physics in our calculations, this 
problem did not persist. 

The most obvious flaw besides the wrong velocity, is the 
steepness of the continuum in the wavelength range from 5200 
to 6200 A. After calculating a grid of models it became clear that 
there is no combination of the model parameters that reproduces 
this slope of the continuum. 

However, it is plausible and to be expected that there is ad- 
ditional extinction present in the host galaxy. Therefore, we cor- 
rect the observation for additional reddening. With an assumed 
extinction of E{B - V) - 0.3, the continuum of the model from 
Fig. |2]was matched very well (See Fig. [Sj. Since the assumed 
value of extinction is an arbitrary choice and is only a lower 
boundary for possible values of the extinction we ran a number 
of tests to place limits on the extinction coefficient. Large values 
of E{B-V) demand high model temperatures, but higher temper- 
atures yield synthetic spectra that are dominated by the Balmer 
series of hydrogen and other spectral features were diminished. 

Hence we decided to choose the extinction as small as 
needed to match the model. During the progress of the analy- 
sis, due to extensive NLTE calculations, non-solar abundances 
and taking the other observed epochs into account. In the final 
model the extinction was chosen to E{B - V) = 0.4. Since the 
continuum wasn't a reliable temperature indicator anymore we 
had to rely on reproduced spectral details to determine the model 
temperatures. 

Because the model spectrum in Fig. |2] was very promising, 
we have based our further modeling on the parameters of this 
model. It already includes hydrogen, helium and the first four 
ionization stages of iron in NLTE. During the course of the mod- 
eling we included all the species in Table [T] in NLTE. Further 
extensive model grids were calculated, varying the model tem- 
perature, the density exponent, time since explosion (important 
for the radioactive decays), and the velocity field. 



Table 1. This table lists all ionization stages that were treated in 
NLTE. 



z 


Ions 


Z 


Ions 


1 


Hi 


14 


Si i-iii 


2 


He I 


15 


Pi-ii 


6 


Ci-ii 


16 


Si-ii 


7 


Ni-ii 


19 


Ki-ii 


8 


Oi-ii 


20 


Ca i-m 


10 


Nei 


25 


Mn i-n 


11 


Na i-ii 


26 


Fe i-iv 


12 


Mg i-iii 


27 


Co i-m 


13 


Ah-ii 


28 


Ni i-iv 



We found that the spectrum of Feb. 4th was best fitted with 
an effective temperature of 5800 K and a density exponent n - 9. 
The exact age of the SN proved to be not particularly important. 
However, it is clear that the radioactivity is necessary to repro- 
duce the features in the observed spectra as it influences the tem- 
perature structure. 

To further improve the model quality we dropped the as- 
sumption of solar abundances and did an abundances analysis. 
Assumed under abundances of specific species such as iron, car- 
bon, or nitrogen of up to [Z/H] = -6 resulted in little or no effects 
on the spectra. The most noticeable effect resulted from an over- 
all under abundance of all the metals. The spectra fit significantly 
better if the metalicity was [M/H]= log ^ ~ -0.5. However, 
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Fig. 2. The plot shows the first model that reproduced some of 
the observed spectral features. No model was found that matched 
the observed rise of the continuum, hence the extinction inside 
the parent galaxy appears to be important. The spectra were nor- 
malized in such a way that they match the continuum above 6900 
A. 
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Fig. 3. The same model as in Fig.|2]is shown. The observation 
was dereddened using a higher extinction of E(B - V) = 0.3. The 
spectra were normalized to match the continuum above 6900 A. 



since the fit is not good enough to use quantitative methods this 
is only determined by eye. Lower metalicities seem to fit the 
observation for a broader range of values. Hence the progenitor 
seems to have been a metal-poor star. 

The variation of the velocity field was required in order for 
the models to match the observations. Gradually changing the 
velocity field resulted in a different wavelength shift and only 
slightly changing spectral features. No particular velocity could 
be favoured other another in a range of a few hundred km s ' . 

The variation of the velocity was problematic in the first 
models due to the fact that slower models didn't converge prop- 
erly or didn't reproduce the spectra as well as before. However, 
due to the refining of the atmosphere structure via inclusion 
of more NLTE species the models became more robust to the 
change of the velocity parameter 

In Figs. |5}{8] we show the final model spectra for all four 
observed epochs. In all Figs, the observations were corrected for 
an extinction of E(B - V) = 0.4. This change of the extinc- 
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Fig. 4. A model is compared to the observation on Feb. 4. Solar 
abundances are assumed. The other key-parameters are Tmodei - 
5800 K and vo = 4900 km s"'. The observation was corrected 
for E{B - y) = 0.4. The spectra were normalized by area. 
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Fig. 5. The best fitting synthetic spectrum is compared to the 
observation on Feb. 4. All model parameters are the same as 
in Fig. |4] with the exception of the metalicity which is [M/H]x! 
-0.5. The observation was corrected for E{B -V) - 0.4. In com- 
parison with the model from Fig. |4] the fit is now much better. 
The spectra were normalized by area. 



tion from the initial guess is mainly due to the use of non-solar 
abundances. The under-abundance of the metals increased the 
flux in the blue part of the spectrum and therefore the observa- 
tions demanded more de-reddening. In order to support the large 
reddening that we find, we have searched for Na ID interstellar 
absorption lines in the host galaxy. Since there is a strong broad 
SN feature just as the wavelength of the Na ID lines they don't 
clearly stand out, but there is a hint of an extra narrow absorption 
line at the right wavelength. Thus, our somewhat high extinction 
value is reasonable. 

The important model parameters of the different observa- 
tions are summarized in Tabled 

5. Conclusion 
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Fig. 6. The best fitting synthetic spectrum for Feb. 5 is com- 
pared to the observation. The model has rmodei = 5700 and 
Vo = 4850km s"' . The spectra were normalized by area. 
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Fig. 7. The best fitting synthetic spectrum for Feb. 9 is com- 
pared to the observation. The model has Tinodei = 5600 and 
v'o = 4500km s"' . The spectra were normalized by area. 
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Fig. 8. The best fitting synthetic spectrum for Feb. 1 1 is com- 
pared to the observation. The model has rmodei = 5400 and 
Vo = 4300km s"' . The spectra were normalized by area. 





We have modeled the early spectra of the Type II supernova SN 
2003Z. The results show that there is substantial extragalactic 
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Table 2. Summary of the basic parameters T„,odei and velocity 
for the different days. There is no error estimation for the ve- 
locity as a range of velocities fit the spectra and the given value 
was just picked to match the observed redshift. The models fit 
best with metalicities [M/H]!^; -0.5, however, there is quite a 
range of possible values that cannot be ruled completely out, so 
[M/H]= -0.5 + / - 0.3/0.5. The luminosities are relatively well 
determined - Armodei - ±100K - for the assumed extinction, 
but depend strongly on this parameter The pressure po is deter- 
mined at T = 1 in the radial r-grid. 



February 




[M/H] 


vo 


n 


Po 




[K] 




[km s-'] 




[10-** bar] 


4th 


5800 


-0.5 


4900 


9 


4.43 


5th 


5700 


-0.5 


4850 


9 


5.41 


9th 


5600 


-0.5 


4500 


9 


6.83 


11th 


5400 


-0.5 


4300 


9 


9.76 



Zel'dovich, Y. B. & Raizer, Y. P. 1967, Physics of shock waves and 
high-temperature hydrodynamic phenomena (New York: Academic Press, 
1966/1967, edited by Hayes, W.D.; Probstein, Ronald F.) 



reddening towards the SN, probably due to dust in the plane of 
its parent galaxy. For each observed epoch we have determined 
the best-fit parameters of the NLTE models as summarized in 
Tabled 

The best fits require sub-solar metalicities of about [M/H]^ 
-0.5, indicating that the p rogenitor of SN 20 03Z was a quite 
metal-poor star (see also iBaron et akl (l2003h ). Unfortunately 
there are only very few observations available for this object and 
in particular, no accurate light curve is known, therefore we have 
assumed that the relative fluxes are accurate. Without photome- 
try we can not calibrate the relative fluxes and thus our results 
are sensitive to the accuracy of the relative flux calibration. 
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